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Abstract. Using a non-perturbative method developed in a previous work (paper II), we derive the probability 
distribution V(Sr) of the density contrast within spherical cells in the quasi-linear regime for some non-Gaussian 
initial conditions. We describe three such models. The first one is a straightforward generalization of the Gaussian 
scenario. It can be seen as a phenomenological description of a density field where the tails of the linear density 
contrast distribution would be of the form Pl(Sl) ~ e~' SL ' , where a is no longer restricted to 2 (as in the 
Gaussian case). We derive exact results for V(Sr) in the quasi-linear limit. The second model is a physically 
motivated isocurvature CDM scenario. Our approach needs to be adapted to this specific case and in order to 
get convenient analytical results we introduce a simple approximation (which is not related to the gravitational 
dynamics but to the initial conditions). Then, we find a good agreement with the available results from numerical 
simulations for the pdf of the linear density contrast for Sl,r <; 0. We can expect a similar accuracy for the 
non-linear pdf V(5r). Finally, the third model corresponds to the small deviations from Gaussianity which arise 
in standard slow-roll inflation. We obtain exact results for the pdf of the density field in the quasi-linear limit, to 
first-order over the primordial deviations from Gaussianity. 
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1. Introduction 

In usual cosmological scenarios, large-scale structures in 
the universe are generated by the growth of small pri- 
mordial density perturbations, through gravitational in- 
stability. At large scales or at early times one can use a 
perturbative approach to obtain the first few moments 
of the evolved density field. Moreover, as described in 
Bernardeau (1992) and Bcrnardeau (1994a) one can ac- 



tually sum up the perturbative series at leading order in 
the limit a — > (where a is the rms density fluctuation) 
to obtain all order cumulants. This yields the probability 
distribution function (pdf) V(5r) of the density contrast 
within spherical cells. However, this method only applies 
to Gaussian initial conditions. Moreover, it may miss some 
non-perturbative effects. 



In a previous paper (paper II) we developed a non- 
perturbative method to derive the pdf V(Sr) in this quasi- 
linear regime. It is based on a steepest-descent approxi- 
mation which yields asymptotically exact results in this 
limit. This allowed us to justify the results obtained by 
perturbative means and to correct some misconceptions 
related to non-perturbative effects. Another advantage of 
our approach is that in principle it can also be applied 
to non-Gaussian primordial density fluctuations. However, 
the determination of the relevant saddle-points needed in 



order to get simple analytic results may then be more dif- 
ficult. In particular, the method may need to be adapted 
to specific cases. 

Thus, in this article we show in details how to ob- 
tain the pdf V(Sr) for three specific non-Gaussian models, 
in the quasi-linear regime. Indeed, although observations 
are consistent with Gaussian initial conditions so far, it 
is of interest to investigate a broader class of models un- 
til any non-Gaussianity is definitely ruled out. Moreover, 
there exist some reasonable physical scenarios (though 
somewhat more contrived than the standard CDM model) 
which give rise to non-Gaussian primordial fluctuations. 
Finally, even standard inflationary scenarios lead to small 
deviations from Gaussianity for the primordial density 
fluctuations. 

This article is organized as follows. First, in Sect .0 we 
recall the path-integral formulation which allows one to 
write an explicit expression for the pdf V(Sn) in terms of 
initial conditions. Then, in Sect]3 we apply our method 
to a non-Gaussian model which is a straightforward gen- 
eralization of the Gaussian case and could be used as a 
phenomenological tool to reproduce observations in case 
some departure from Gaussianity would be measured (e.g., 
in the CMB data). Then, in Sect. |] we investigate a phys- 
ically motivated model which describes an isocurvature 
cold dark matter scenario ( Peebles (1999a) ). In particu- 
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lar, we obtain a good agreement with the available re- 
sults from numerical simulations for the pdf of the lin- 
early evolved density field. Finally, in SectJq we consider 
the small primordial deviations from Gaussianity which 
arise in standard slow-roll inflation. 

2. Generating functions. Gaussian initial 
conditions 

In this article, we investigate the probability distribution 
function (pdf) V(5 R ) of the density contrast S R within a 
spherical cell of comoving radius R, volume V: 



6r = 



d 3 x 

V 



<J(x) 



(1) 



Here 6(x) is the non-linear density contrast at the comov- 
ing coordinate x, at the time of i nterest. I n this section, 
following the method developed in paper II , we recall how 
one can express the pdf V(S R ) as a function of the initial 
conditions through a path-integral formalism. In order to 
introduce our approach we first briefly consider the case of 
Gaussian initial conditions which is most familiar. In the 
next section, we extend our formalism to a non-Gaussian 
model which can be seen as the simplest generalization of 
the usual Gaussian case. 

Rather than trying to directly evaluate the pdf itself, it 
is actually more convenient to study its Laplace transform 
i/j(y) given by: 



i>(v) 



,-vSr\ 



dS R e-* s * V{5 R ). 



(2) 



Here, the symbol (..) expresses the average over the initial 
conditions. Then, the last term in eq.(g) can also be seen 
as the definition of the pdf V(8r). Indeed, the pdf can be 
recovered from ip(y) through the standard inverse Laplace 
transform: 



V(S! 



-\-ioo 



dy_ 
2iri 



?ySa 



ip(y)- 



(3) 



Moreover, the generating function i/j(y) also yields the mo- 
ments (S R ) through: 



i>(y) = 



E 

9=0 



{-y) q 



(s q R : 



(4) 



Then, we need to compute ip(y) as an average over the 
initial conditions, using the first equality in eq.(0). To do 
so, one simply needs two pieces of information. First, we 
must know the functional <5r[<5l(x)] which yields the exact 
non-linear density contrast 5 R over the cell V which arises 
from the gravitational dynamics of the linear density field 
(5l(x). Indeed, as shown for instance in paper I the initial 
conditions can be defined by the linear density contrast 
5l(x) where we only keep the linear growing mode. This 
does not assume that the exact non-linear density field 
(5(x) can be written as a series expansion over powers of 
the linear fi eld <$r,(x ) . In fact, such a series is only asymp- 
totic (e.g., paper I , paper V ). Hereafter, we note S L (x) 



the growing mode of the linear density contrast at the 
time of interest, where we compute V{5 R ) (i.e., to sim- 
plify notations we do not write explicitly the time depen- 
dence). Second, in order to perform the average in eq.pl) 
we need the weight which is associated with all possible 
fields <5l(x). In other words, we must specify the probabil- 
ity distribution of the random field <5l(x). Thus, the first 
point expresses the physics of gravitational interactions 
while the second point describes the initial conditions of 
the system. 

In the case of Gaussian initial conditions, the statistics 
of the random field 8l (x) are fully defined by the two-point 
correlation: 



A L (xi,x 2 ) = (<5 L (xi)5 L (x 2 )). 



(5) 



The kernel A^, is symmetric, homogeneous and isotropic: 
Ai(xi,x 2 ) = Al(|xi— x 2 |). It is convenient to express A^ 
in terms of the power-spectrum P(k) of the linear density 
fluctuations. To do so, we define the Fourier transform of 
the density field as: 




dke 



■ik.x 



*(k) 



(2tt) s 



<f(x) 



(6) 



Next, we define the Fourier transform of the kernel A^ by 
the property: 



/: 



.A L ./ 2 = y"dkxdk 2 / 1 (k 1 )*.A L (k 1 ,k 2 )./ 2 (k 2 ) 



(7) 



for any real fields j\ and / 2 , where we introduced the 
short-hand notation: 

/i-A L ./ 2 = /dxxdx, / 1 (x 1 ).A L (x 1 ,x 2 )./ 2 (x 2 ). (8) 

Using eq.(g) this implies: 

A L (k 1; k 2 ) = /dxxdx, e *(k2.x 2 -k 1 .x 1 ) Ai(xi,x 2 ) (9) 



which gives: 

A L (k 1; k 2 ) = (2nf P(h) Mki - k 2 ) 



(10) 



where we defined the power-spectrum P(k) of the linear 
density contrast by: 

((J i (k 1 )<J L (k 2 )) = P(fci) Mki + k 2 ). (11) 

Finally, the inverse A^ 1 of the kernel A^ is (see paper II): 



AZ 1 (k 1 ,k 2 ) = 



S D (ki -k 2 ) 



(12) 



P(ki 

which implies that A^ is positive definite since we have 

IMk)| 2 



6 L .Al\S L = I dk 



P(k) 



(13) 



where we used 8l{— k) = <5j,(k)* for real fields <5l(x), see 
eq-i). 
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Thus, for Gaussian initial conditions the statistical 
properties of the random field <5l(x) are determined by 
the kernel Al given in eq.(lQ) (or equivalently by A^ 1 
given in eq.(n2|)). In particular, the average over the initial 
conditions in cq.(pl) can be written as the path- integral: 



^( 2 /) = (DctA^ 1 ) 1/2 |[d ( 5 L (x)] 



-vSr[S l 



\Sl.A; 



,(14) 



This expression merely means that in order to compute 



the average (e 



-v$r\ 



we simply need to sum up the contri- 



butions of all possible linear density fields 5l (x) to which 
we associate a Gaussian weight which is proportional to 
e ~(8 L .A- 1 .s L )/2_ The norma ii zat i on f actor (DctA^ 1 ) 1 / 2 en- 
sures that V'(O) = 1; as implied by the definition (0). 
Here DetA^ 1 is the determinant of the kernel A^ 1 . Thus, 



eq.(14) yields an explicit expression for the Laplace trans- 
form ij)(y) (though it may be difficult to obtain numer- 
ical results from such a path- integral). This provides in 
turn the pdf V{5r) through the inverse transform (0). In 
paper II we showed how to get the generating function 
ip(y) from eq.dl4) in the limit -P(fc) — ► (i.e. for small rms 
density fluctuations) using a steepest-descent method. In 
the next section, we apply this method to a closely related 
non-Gaussian model. 

3. Generalization of the Gaussian weight 

In this section, we describe a non-Gaussian model which 
can be seen as the simplest generalization of the Gaussian 
scenario. This model is not derived from physical prin- 
ciples. It simply provides a phenomenological tool which 
can describe initial conditions such that the tails of the 
linear pdf Vl(8l) of the density field are of the form 
Vl(Sl) ~ e~' SL \ , where a is no longer restricted to 2 (as 
in the Gaussian case) . In addition, it allows us to show on 
a s imple ex ample the power of the method we developed 
in [paper II and how it can be applied to non-Gaussian 
initial conditions. 



3.1. A non-Gaussian model 

In the case of Gaussian primordial density fluctuations, 
the average involved in the definition (H) of the generating 
function ip(y) (i.e. the Laplace transform of the pdf V{8r)) 
is simply given by the Gaussian weight exp(— |#£,.A£ .#£,) 
in the path-integral (S3) which takes the average of e~ v&R 
over all possible initial states, as we recalled in the pre- 
vious section. A straightforward generalization of the 
Gaussian case (^4|) is to modify this weight. Thus, one 
can consider the case where the average over the field 5l 
is given by: 



?/%) = M / [(W L (x)] e - y5R ' w ^ 



with: 
W[S L ]=J2*iV L A£-8L) ai/2 



;v 



(15) 



(16) 



where ot\ < .. < a^ and A at > 0. Here TV is a normal- 
ization constant, so that ip(0) = 1, and we still take the 
kernel A^, 1 to be of the form (O) hence A^, 1 is positive 
definite. We also define a quantity aw in terms of Ayp by: 



, ff (ijfE0^A ff ( Xl ,x 2 ). 



(17) 



However, contrary to the Gaussian case we now have: 
(6l(xi)6Ux-2)) ± A w ( Xl ,x 2 ) , (Sl >R ) ? a w (R) 2 . (18) 

The form ( J15| ) corresponds to a large density contrast tail 
of the form Vl(Sl,r) ~ e _ '* x,,B ' N for the linear density 
field. 

In order to compute the Laplace transform ip{y) we 
simply follow the steps of the calculation developed in 



paper II. Thus, we first define a new generating function 



tp{y) through: 



4>(y) 



-ip(ya 



,cr w )/a v 



(19) 



in order to factorize the amplitude of the kernel A\y. This 
yields from eq. (|l5|) : 



-f(v,"w)/^ = tf I [dS L (x)] e~ Sw/ 



where the "action" 5V[<y is given by: 



A 



S w [Sl] = y S R [8 L ] + J2 x > ^(Sl-A^.S^ 



i/2 



(20) 



(21) 



The quasi-linear regime corresponds to the limit aw ~* 
at fixed y. Then, we see that all terms with 1 < i < N — 1 
in the action Sw in eq.(pi|) vanish as ex ct^ - " 1 since 
A^ 1 oc cr^ 2 , see eq.(|l7|). Therefore, in the limit aw — > 
we are left with the action S a : 



S a [S L ] = y S R [S L ] + l^ (S L .A w l .S L ) 



nq/2 



(22) 



1=1 



with a = oln and we normalized Aw such that Ajv = 1/a. 
The fact that the terms with 1 < i < N — 1 disappear in 
the quasi-linear limit expresses the fact that this limit is 
actually a "rare-event limit". Indeed, as seen in paper II 
for the Gaussian case (and this remains valid here) a finite 
y corresponds to a finite density contrast. Then, in the 
limit aw — ► these finite density contrasts (even though 
small) become very rare events. On the other hand, the 
tails of the functional distribution V[Sl] defined by the 
weight VF[^l] are governed by the highest power in eq.(|l6|). 
This implies that in the quasi-linear limit the only relevant 
term is the highest power ajy. Note that we can expect 
the validity of the quasi-linear limit to extend to larger 
values of aw in the case N = 1 than for N > 1. Indeed, in 
this latter case we can expect the terms l<i<iV— lto 
make a significant contribution for aw ~ 1 which is not 
taken into account at all by the quasi-linear limit. 
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3.2. Steepest-descent method 

The action S a is indep endent o f the normalization of the 
kernel A\y- Thus, as in paper II we can apply the steepest- 
descent method in the limit aw — > 0. Indeed, it is clear 
than in this limit the path-integral ( |2C| ) is dominated by 
the global minimum of the action S* Q [(5i] while the contri- 
butions from other points Sl are exponentially damped. 
Then, the steepest-descent method yields asymptotically 
exact results in this limit. We briefly describe below the 
main steps of this steepest-descent method. Since the 



derivation is very close to the one performed in paper II 



for the Gaussian case we refer the reader to that paper for 
the details of the calculation. Note that the Gaussian case 
actually corresponds to a = 2. 

In order to apply the steepest-descent approximation 
we first need to find the global minimum of the action 
SoJiSl]. The condition which expresses that the point 5l 
is an extremum (or a saddle-point) is: 



8S a 



8(8l{x)) 



= for all x, 



(23) 



where <!>/<5(<5l(x)) is the functional derivative with respect 
to 8l at the point x. This constraint also writes: 



fc(xO = -f (fc.A^.fc) 1 -"/ 2 



'w 



/dx" A M ,(x',x" 



S(Sr) 



(24) 



Thus, we obtain exactly the same profile for the saddle- 
point as in the Gaussian case. Indeed, we can see that 
a =/= 2 only modifies a numerical multiplicative factor 
in the r.h.s. of eq.(|24|) (i.e. independent of x' and x") 
which could formally be absorbed into y through y — ► 
ya w ~ a (8L-A w V .£l) 1_Q!//2 . Then, as in paper II we obtain 
a spherically symmetric saddle-point. Note that the ex- 
istence of a spherical saddle-point could be expected a 
priori, since the initial conditions are homogeneous and 
isotropic and we study the density contrast within spher- 
ical cells. Therefore, the very problem we investigate is 
spherically symmetric. However, a priori this spherical 
saddle-point is not necessarily the global minimum of the 
action (in full generality it might as well be a maximum) . 
This will need to be checked a posteriori. 

The fact that we obtain a spherical saddle-point 
greatly simplifies the problem since it means that this 
point Sl is described by the well-known spherical collapse 
solution of the gravitational dynamics. This reads: 



R 



F[Sl,r l } 

-- (1 + 5 R )R 3 



(25) 



where the function T[8l.r l ] is given by the usual spher- 
ical collapse s o lution of the equations of motion (e.g., 
paper II ). 



Peebles (1980) 



The second equation in fl25|) 
merely expresses the conservation of mass and the fact 
that the matter enclosed within the radius R in the actual 
non-linear density field actually comes from a Lagrangian 



comoving radius Rl, for a spherical initial condition. The 
eq.(Eq) provides the functional <5r[<$l(x)] for spherical 
states <5z,(x). As shown in paper II this is sufficient to de- 
rive the spherical saddle-point which obeys the condition 
(p4[). Using the results of paper II with the substitution 



V —> V a w a (^l -Ar/ -Sl) 1 "^ we eventually obtain: 



1L,R L 



(^.A^)(«-2)/2 



'W 



F'[8l,r l ]o w (Rl)/o w (R) 



1 - F [8l.r l 
together with: 



3R 2 r 



1L,R L 



1 dq-yy 

o-w(Rl) dR 



(Rl) 



TL,R> 



J L,R L 



A W (R',R L ) 
A W (R L ,R L )' 



(26) 



(27) 



These two equations fully define the spherically symmetric 
saddle-point (5l(x). The implicit equation ( p6| ) determines 
8l,r l while the radial profile of this initial state is given by 
eq. p7l) . In eq. pTJ) we introduced the spherically averaged 
kernel Ayy (R 1 ,R 2 ) defined by: 



A,, (/(•,. /?_,) = / ^ / ^A ff ( Xl ,x 2 ). 



Vi 



Vi 



V-2 



Vi 



(28) 



As in paper II these equations can be simplified by intro- 



ducing the functions t(8l) and G{t) defined by: 
- 8 L a w (R) 



t(8l) = 



and: 



<j w [{\ + T[8 L ]y/*R} 



G(t)=T[8 l (t) 



(29) 



(30) 



where <5l(t) is defined by eq.(|2^). Using eq.(p9[), we note 
that for a power-law "power-spectrum" P(k) oc k n we 
have aw{R) <x R~( n+3S >/ 2 so that eq.(30) simplifies to: 



Q{j)=T -t{1 + Q\t\) 



-(n+3)/6 



(31) 



Note however that in this non-Gaussian case the "power- 
spectrum" P(k) does not obey eq.(O). Nevertheless, it 
still provides a measure of the power associated with the 
wavenumber k. 

Finally, the limiting generating function ip(y) defined 
by <fi(y) = tpiViV — 0) is given by the value of the action 
S a at this minimum which yields the implicit system: 



IQ-2 



■vQ'(t) 



(32) 



<p(v) = V G( T ) + 



Of course, this result is similar to the Gaussian case which 
can be recovered by taking a = 2. Note that eq.(p2) has 
only been derived for real y (hence r is real), so that |r| 
is the absolute value of r. 

Thus, the steepest-descent method allows us to com- 
pute the Laplace transform ip(y) in the quasi-linear regime 
through eq.(p2|) which yields ip(y) as in cq.(]19|). This will 
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give the pdf V(5r) through eq.(0). However, we first need 
to check that the spherical saddle-point defined by eq. (Eq) 
and eq. (M) is the global minimum of the action S a . For 
small positive y this can be rigorously proved as for the 
Gaussian case, using the fact that the kernel A^ is posi- 
tive definite. We shall not repeat this discussion here (see 



Sect. 3. 4 in paper II). For negative y this cannot be proved 
as it may actually happen that the spherical saddle-point 
is only a local minimum. This also occurs in the Gaussian 
case for a power-spectrum with n < 0. This case is dis- 
cussed in great details in paper II where we show that 



the steepest-descent method remains useful but requires 
a careful justification. Since the same discussion can be 
applied to the non -Gaussian model we investigate here we 
refer the reader to paper II for a description of such cases 
and we shall only give a brief comment on this point in 
the next section. 

3.3. Geometrical construction 



As for the Gaussian case studied in paper II we can give a 
geometrical construction of the generating function (f(y) 
defined by eq.(J32|). Indeed, this latter expression can also 
be written: 



v(y) 



y G(t) + L - L - 
a 



(33) 



The minimum which appears in this expression merely 
expresses the fact that by sheer definition of the steepest- 
descent method the generating function f(y) is governed 
by the global minimum of the action S a . 

Then, we see that the geometrical construction dis- 



played in Fig. 4 in paper II still applies to this non- 



Gaussian case, with the modification that the parabola 
h — r 2 /(2y) are replaced by the curves h— \r\ a /(ay). That 
is, the minimum point r which yields <p{y) in eq. (p3|) is 
given by the first contact of the curve Q{t) with the curves 
h— \r\ a /(ay) of varying height h. For y > we start from 
below at h = — oo while for y < we start from above at 
h = +oo. 

As in the Gaussian case, we find that for some values 
of n the action S a has no global minimum. This can be 
seen from the geometrical construction shown in Fig. 4 or 
simply from cq.(|33|). Indeed, as noticed in |papcr II from 
eq.(|3l|) we have the following asymptotic behaviour for 
Q{r): 



high densities 



-oo, g 



g ~ (-r) 



__>|6/(n+3) 



(34) 



Using this expression we see that for y < the term 
in the brackets in the r.h.s. of eq.(B3) is not bounded 
from below if a < 6/(3 + n), since in this case it goes 
to — oo for r — ► — oo. For the Gaussian case this occurs 



for n < 0. However, as explained in paper II the steepest- 



descent method developed in Sect. 3.2 remains useful and 



the pdf T(5r) is still governed by the spherical saddle- 
point defined by eq.(|26|) and eq. (g7j) . This can be seen 
by the following remark. If we consider the quantity dJ q 



rather than <5r, where q is an odd integer, we can again 
apply the steepest-descent method which yields the asso- 
ciated generating function tp q (y) as in eq. (J33J) . However, 
the function Q q {r) is now given by g q (r) — ^(r) 1 / 9 . Then, 
if we choose a sufficiently large value for q the behaviour 
of the r.h.s. bracket in eq.(|33|) is dominated by the term 
\r\ a at large \t\ so that the action now exhibits a global 
minimum. Of course, this point is still the spherical saddle- 



point we obtained in Sect, 3.2 (since the initial conditions 
5l which give rise to a fixed value of 5r do not depend on 
whether we study 5r itself or S R !). Then the pdf V(Sr) 
could be recovered from V{5^ q ) through a mere change 
of variable. This implies in turn that the pdf V(8r) is 
still governed by the standard spherical saddle-point, even 
though the latter is no longer the global minimum of the 
action. However, as described in paper II the treatment of 



this case requires some additional care and it is actually 
associated with a break-up of perturbative theory. 

3.4. Calculation of the pdfV{5 R ) 

Finally, we describe how to compute the pdf V{8r) itself 
from the results obtained in Sect. 3.2. We consider the case 
n = — 1 and a = 3 where the action S a always shows a 
global minimum. Then, the Laplace transform ifj(y) is fully 
defined by the generating function tp(y) given by eq. (|32|) . 
In order to compute the pdf from eq.(p3) we simply 



use the inverse relation (0) which yields: 



V(S R ) 



-4-zoo 



dy 



Av s r-v(v)]/cw 



-loo Altl/Cyy 



(35) 



Then, to perform the integration over y we need the ana- 
lytic continuation of tp(y) over the complex plane. Indeed, 
note that eq.(|32]) and eq.(J33) were actually derived for real 
y (hence t was real). In particular, we must continue the 
absolute value \t\ over the complex plane. Moreover, we 
need to specify the integration path over y (since ip(y) will 
usually be singular at the origin). First, the integration 
path starts from the real axis at the saddle-point (r c ,y c ) 
given by: 



(j/ c ) = o with x(y) = Srv - tp(y)- 



dx 
dy 



From cq.(|32j) we have: 

G(i~c) = Sr since <p'(y) =G(t). 



(36) 



(37) 



Thus, as for the Gaussian case, we see from eq.(30) that 
the triplet (Sn,T c ,y c ) is also the triplet (8r, t, y) we ob- 
tained in Sect. [3.2| to get tp(y). This simply means that 
V(6it) at the point Sn is governed by the neighbourhood 
of the saddle-point <5l(x) obtained in Sect.^, which obeys 
J~{^l,Rl) — Sr. This is actually quite natural. Then, we 
need to continue the function (p(y) over the complex plane 
from the neighbourhood of the point y c on the real axis. 
Thus, if t c > (i.e. y c > 0) we replace \t\ in eq.(32) by r, 
while if t c < (i.e. y c < 0) we replace \t\ by (— r). Then, 
we have the usual analytic continuation of the power-laws 



(i 
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r a or (— r) Q . The argument of these power-laws is real 
positive at the starting point (r c ,y c ). This provides the 
analytic continuation of the generating function <p(y). The 
function Q(t) is also analytic around the real axis because 
the spherical collapse solution F(8l) can be expressed in 
terms of analytic functions (e.g., trigonometric or hyper- 
bolic functions for a critical density universe). 

Next, the integration path in the complex plane fol- 
lows the steepest-descent contour which runs through the 
saddle-point y c . It is given by the constraint Im(x) = 
and it is orthogonal to the real axis at the point y c . This 
path is also symmetric about the real axis which implies 
that the result for V(8r) is real. Note that in the quasi- 
linear limit aw —* the contributions to the integral 
( p5| ) only come from an infinitesimal neighbourhood of 
the saddle-point (r c ,y c ) around the real axis. 

Usually the function tp(y) is not regular at the origin 
since we have <p(y) — — (1 — l/a)\y\ a ' i - a ~ 1 ' > for small real y. 
Then, the pdf V(8r) is not analytic at Sr — and the mo- 
ments of the pdf cannot be recovered from the expansion 
(Q). This is not the case for Gaussian initial conditions 
where both ip(y) and V(Sr) are regular at the origin. 




Fig. 1. The pdf V(Sr) for n = -1, fl m = 1, a = 0.74 
and a = 3. This corresponds to aw — 0.65. The solid line 
shows the theoretical prediction from (|32|). The dashed- 
line displays the linear pdf Vl(8l,r) of the linearly evolved 
density field. The dotted-lines show the linear and non- 
linear pdfs obtained for Gaussian initial conditions with 
the same linear rms density fluctuation a (the non-linear 
pdf is the asymmetric curve with the extended high- 
density tail). 



As an example, we show in Fig.n] the pdf V(5r) for 
the case a = 3 and n = — 1 (solid line). We also have 
a = 0.74 where a 2 = (8\ R ) is the rms density fluctua- 
tion. Since a ^ aw we numerically compute a from the 
second moment of the linear pdf Vl ■ Indeed, the steepest- 
descent method developed in Sect. 3^ also yields the linear 
pdf Vl{Sl,r) °f the linearly-evolved density field. In this 
case we have Sr = 5l and R^ = R, hence we simply need 
to take Q{t) = — r in eq.(p2[) (of course, for a = 2 this 



gives back the Gaussian) . The non- linear evolution of the 
density field increases the value of ^(Sr) at large posi- 
tive Sr (it actually spreads the linear pdf towards larger 
value of the density) as in the Gaussian case but the cut- 
off remains much sharper. It also induces a cutoff at low 
densities Sr> —1 which expresses the fact that the actual 
density p is positive. 

For comparison, we also display in Figffl the results we 
obtain for Gaussian initial conditions with the same linear 
variance a (dotted lines). These curves were derived in 
paper II. As could be anticipated from eq.(15) the large 
density contrast cutoff of the linear pdf is much sharper 
than for the Gaussian case: it goes as V(Sl,r) ~ e> SL > H > 
since here we take a = 3. Moreover, we can check that 
for the actual non-linear pdfs this trend remains valid. 



More precisely, as shown in 3aper II the high-density tail 



of the non-linear pdf for Gaussian initial conditions is of 
the form: 



Gaussian scenario : 5r 3> 1 : ^(Sr) 



-8 



(n + 3)/3 



1^ 



(38) 



while for non-Gaussian initial conditions we obtain: 



non-Gaussian : 6r > 1 : V(S R ) - e^^ 3 ' 76 /^. (39) 



These behaviours are obtained from eq. 
eq.(|37|) which yield: 



V(S R ) 



eq.(|35j) and 

(40) 



while the asymptotic behaviour of t c (Sr) for large Sr is 
given by eq.(^J) and eq.(p7|). 

As noticed above, the non-Gaussian model ( |l6| ) is not 
derived from a physical scenario of the primordial uni- 
verse. It should be viewed as a simple model to describe 
initial density fluctuations which exhibit a non-Gaussian 
tail. Then, the analysis performed in the previous sections 
shows how one may derive in a rigorous manner the grav- 
itational dynamics of density fluctuations in the quasi- 
linear regime. Moreover, we can expect the form (39) of 
the high-density tail to apply to any model which obeys 
~P(^l,r) ~ e~' l5l " R ' at large densities, whatever the de- 
tails of the model. In fact, we have shown that in the 
quasi-linear regime the tail (p59) (actually the generating 
function (p(y) itself !) is common to all models of the class 
described by eq.(|lq) with o.n = a. 

In practice, since observations have not shown any de- 
viations from Gaussianity so far, the parameter a is con- 
strained to be close to 2 (the Gaussian value). Note that 
this family of models presents the advantage of containing 
the usual Gaussian scenario as a particular case. 

4. An isocurvature CDM model 

In addition to generalized non-Gaussian models as the one 
studied in Sect. g which are rather ad-hoc tools, there ex- 
ist some specific non-Gaussian models which arise from 
a physically motivated description of the primordial uni- 
verse. These scenarios are of great interest as they could 
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provide a reasonable alternative to the standard simple in- 
flationary model. Moreover, it is prudent not to disregard 
a priori sensible models of structure formation in view of 
the rather indirect character of the probes of the early 
universe which are available to us. 

4.1. Initial conditions 

Thus, in this section we consider a specific isocurva- 
turc cold dark matter model which was presented in 
Peebles (1999a) . This is an inflationary scenario which 
involves three scalar fields and it gives birth to non- 
Gaussian isocurvature CDM fluctuations. We refer the 



reader to Peebles (1999a) for a description of the physi- 
cal processes which give rise to this scenario. Then, at the 
end of inflation (at time L) the CDM mass distribution is 



(Peebles (1999a)): 



pi(x) oc ^(x) 2 



(41) 



where <fo(x) is a Gaussian random field with zero mean 
and power-spectrum P ( p i : 

(&(ki)&(k a )) = P 01 (fci)fe(k 1 +k 2 ). (42) 

On the scales of interest for structure formation we have: 



P^(k)cxk n * , n^ 



-2.4 



while below the coherence length x c (i.e. for k 3> l/x c ) the 
power-spectrum decreases as P^k) oc k s with s < —3. 
This coherence length ensures that (<fri(x) 2 } is finite, as 
required by eg. (El] ) . In the case of the model described in 
Peebles (1999a)| we have x c ~ 10 pc. In the linear regime 
below the horizon the CDM density perturbations grow as 
<5j,(x, t) oc D + (t) (we only keep the growing mode) hence 
we write: 



6 L (x,t) = D + (t) (4>i(x) 2 - l) 



(44) 



where we normalized <fo(x) by (</>i(x) ) = 1 and the grow- 
ing mode by D + (ti) = 1. Note that in eg.(B3) and in the 
following we only consider "small" comoving scales x ^ 14 
Mpc, where t he isocurvature CDM transfer function is 
close to unity ( Peebles (1999b) ). On larger scales the den- 
sity fluctuations are damped with respect to eg.(]44|) and 
we should take into account the decrease of the trans- 
fer function. Since we are interested in the statistics of 
the density field at an arbitrary time t it is convenient 
to rescale the Gaussian field <\>i. Thus, we define </>(x) = 
£> + (£) 1/2 <^(x) so that eg.(E|) now reads: 



<5 L (x) = 0(x) 2 



at the redshift of interest. The power-spectrum P$(k) is 
still of the form (M3f). Note that the relation between 
the fields <5i(x, t) and (ft(x, t) depends explicitly on time 
through the last term in eg. (|45[ ) (to simplify notations we 
usually do not write explicitly the time dependence of the 
various random fields). Indeed, we must recall here that 
(<j) 2 ) is not the spatial average of <fi(x) 2 but the mean over 



all realizations of the Gaussian random field <f>. Therefore, 
it is independent of the peculiar realization <f> which gives 
rise to a given linear density field through eg.pq) (even 
though by ergodicity both averages happen to be egual). 
We define the two-point correlation of the field </>(x) by: 

A (x 1 ,x 2 ) = ( ( /)(x 1 )0(x 2 )) (46) 

and the linear variance over the scale R by: 

a 2 (i?)^^^A,( Xl ,x 2 ). (47) 

Note that eg.fiq) also implies that we can define the av- 

2 > by: 



erage ( 

/a2\ — 



A^(x,x) 



(48) 



which does not depend on x since the initial conditions 
are invariant through translations. From eq.fiq) and the 
fact that </>(x) is Gaussian one obtains the variance cr(R) 
of the linearly-evolved density contrast as: 



a 2 (R) = (8l R ) = 2 



dxi dx 2 



A (xi,x 2 )' 



(49) 



since the two-point correlation A^ of the linear density 
field is given by: 



(43) A L (x 1; x 2 ) = ( ( J L (x 1 )J L (x 2 )) - 2A ( Xl ,x 2 ) 



(50) 



From eg.(J47|) and eg.(|g) we see that for -3 < n^ < -3/2 
we get a(R) oc <r?(i?) and the slope n of the linear power- 
spectrum P(k) of the density fluctuations is related to n^ 
by: n = 3+2ri . This implies n^ ~ —2.4 as in eg. (|43|) in or - 
der to match observational constraints (Peebles (1999b)). 
For numerical calculations we shall use: 



n^ = —2.4 and n = 3 + 2n = — 1. 



4.2. Laplace transform 



(51) 



We again define the generating function ip(y) as in eq. (j2[) 
but the Gaussian average over 6l(x) which appeared in 
eg.(|lj) must now be replaced by a Gaussian average over 
<t>(x): 



My) =Af /[d<Mx)] e -AW4fi/-*. 



(52) 



Thus, our goal is now to estimate the Laplace transform 
?/>(y). To do so, we shall not introduce a rescaled generat- 
ing function <p(y) as in eg. (03) and we work directly with 
ip(y). First, using eg. (45) it is convenient to introduce the 

through 



(45) linear density field 5l(x) into the expression 



the Dirac functional 8d[{4> ) + ^l(x) — </>(x) z 



tl>(y) - / [d0(x)][d<5 L (x)] 6 D [(<f) + 5 L (x) - 4>{*Y 

-yS R lS L ]-±4,.A-\<f> 



x e 



(53) 



Here and in the following we do not write the normaliza- 
tion constant of the path-integral (it will be recovered in 
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the final expressions). Then, using the integral represen- 
tation of the Dirac functional we obtain: 



1>(y) = /[#(x)][dMx)][dA(x)] 

-y6 R [S L ]+i\.{('P 2 )+5L-<p 2 )-i<p.A- 1 .4, 



x e 



(54) 



where A(x) is an auxiliary real field. Note that the expo- 
nent is actually quadratic over </> and we can write eq. ( p4[ ) 
as: 

1>(y) = J[dct>(x)}[d6 L (x)][dX(x)} 

x e -y6 R [6 L ]+i\.((4, 2 )+S L )-^.(A- 1 +2iA).<p / gg s 

where we introduced the kernel: 
A(xi,x 2 ) = A(xi)(5 D (xi -x 2 ). 



Then, the Gaussian integration over 
ward and it yields: 

</%) = |[d5 L (x)][dA(x)][Dct(A^H 
x e -y<5 R -MA.«0 2 )+<5 L )^ 



(56) 
(x) is straightfor- 



■2*A) 



-1/2 



(57) 



Note indeed that the factor (0 2 ) in eq.(p5|) does not de- 
pend on the random field </>(x) over which we integrate, as 
shown by eq.([48|). Next, multiplying the r.h.s. of eq.([57|) 
by (DetA^) -1 '^ (i.e. we change the implicit normaliza- 
tion constant of the path-integral) and using the relation 
ln(DetM) = Tr(lnM) we get: 



<■(.,/) = / [dfo(x)][dA(x)] e -i^MWiAA,) 

-yS R +i\.((<t> 2 )+S L ) 



x e 



(58) 



At this point, it is instructive to consider the linear 
regime when the density field is linearly evolved: <5(x, t) = 
<5l(x, t) oc D + (t). In this case, we have: 



yS R = y 



e(x < R) 

dx S L (x) 



(59) 



where 9{x < R) is a top-hat with obvious notations. Then, 
the integration over <5_l(x) is straightforward and it yields 
the Dirac functional <5d[A(x) + iyO(x < R)/V\. The inte- 
gration over A(x) is now obvious and we obtain for the 
generating function ipL(y) which describes the linear den- 
sity field (again, the subscript "L" refers to "linear"): 



^l{v) 



o-v&l 



o y(0 2 )-iTrln(l+2iA y A ) 



where the matrix A y is given by: 



A y (xi,x 2 ) 



-ly 



(xi < R) 
V 



fo(xi - x 2 ). 



(60) 



(61) 



The normalization of the result derived from the path- 
integral in eq.(f30|) is obtained from the condition V'l(O) = 
0, as implied by the definition oiipL(y) (since (1) = 1). As 
shown in eq.(Q) the generating function ipL(y) gives the 
moments (S q L R ) of the pdf Vl(Sl,r) of the linear density 



contrast. However, it is usually convenient to introduce 
the generating function < &l(2/) of the cumulants (Sl R ) c 
which is given by (see any textbook on statistical theory) : 



Mv) = W^(i/)] = E4^,h> 



9=1 



Thus, from eq. (p0[) we obtain the exact expression: 

My) = (<t> 2 )v - \ Tr Mi + 2 ? A y A ). 



(62) 



(63) 



One can easily check by expanding the logarithm over y 
that $l(j/) is actually quadratic over y (see also eq.([75|) 
below). Hence (5l,r) — as it should. Moreover, we can 
check that we recover eq. (|49| ) for the variance of the linear 
density contrast. 

4.3. Saddle-point contributions 

In the non-linear case it is not so easy to compute the 
expression ( p8[ ) since the integration over 5l (x) no longer 
yields a Dirac functional. Besides, even for the linear den- 
sity field the expression (p3^ is not very convenient since it 
is not obvious to compute the logarithm. Thus, for prac- 
tical purposes it only provides the first few cumulants of 
Sl,r, which can be obtained by expanding the logarithm 
as a series over y and computing the trace of each power of 
A y A,p. Hence we can try a steepest-descent approximation 
to the path- integral (pq), in the spirit of the calculations 
presented in paper H for the Gaussian case or in Sect.Hfor 
a simple non-Gaussian model. More precisely, we may try 
to evaluate the contributions of all possible saddle-points 
to the path- integral (p8|). 

First, we can look for the saddle-points of the exponent 
with respect to <5l(x). This yields: 



iA(x) = y 



S(S R ) 
S(6 L (x)Y 



(64) 



Since the Gaussian random field </>(x) is homogeneous and 
isotropic we still have a spherical symmetry and we can 
again look for a spherical saddle-point. Then, taking the 
mean of eq.(|64|) over a spherical cell V we get: 



, dx' 5(S R ) 
lX *'=vL .177 



v , V S(S L {yf)) 



(65) 



From the very definition of derivatives we can also write 
cq.(J65|) as: 

i\ R , = lim - {S R [6l,r» + eA x (R', R")] - 5 R [6 L , R »]} 

(66) 

where R" is a dummy variable and we introduced the ker- 
nel: 



A^jO.jf £!£<*> .J**"" ■''■ 



^ifR" >R' 



(67) 
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We again define the linear density contrast Sr l and the 
Lagrangian mass scale Rl as in eq. (Eq) and we obtain: 



i>^R' = y 



r[6 L , RL ]A x (Bf,R L ) 



1-F'16l,r l ] 



R 3 ^5 

L,R L \m — 



dR" 



(68) 



in a fashion similar to the calculation presented in paper II 



for the Gaussian case. The only term in the r.h.s. of eq. (|68[) 
which depends on R' is the factor A\(R',Rl) hence we 
have: 



\ R ,ocA x (R',R L ) 

which yields the profile: 
9{x < Rl) 



A(x) = A 



V L 



(69) 



(70) 



which also defines the normalization A . Note that this 
profile for the field A(x) is of the same form as the one 
obtained for the linear regime which involved the Dirac 
functional 5l>[X(x) + iy9(x < R)/V], see the derivation of 
eq-©. 

Then, for such a state A(x) we can estimate the trace 
which appeared in eq.(p8|) as follows. Let us note this trace 
as T(Aq). Then, we have: 



T(A ) = Tr ln(l + 2zAA ) 
with: 



~ (-1)9+1 
9=1 



<1 



(71) 



T„ 



Tr (AA*)« 
/ dxi-.dx, (AA )(xi,x 2 )..(AA )(x g ,xi) 



(72) 



using the definition of the trace: TrM = f dxil/(x,x). 
Next, from eq.(56) we see that for the spherical state ( |70|) 
the matrix (AA^) is given by: 

(AA )(x 1: x 2 ) = A 0{Xl < Rl) A ( Xl ,x 2 ) (73) 

VL 

which yields: 

Tq = K Sv if -^^^I'^-^^^i)- ( 74 ) 

The first term T\ (i.e. q = 1) is simply: 

Ti = A / ^A4x 1 ,x 1 ) = A o <0 2 ) (75) 

where we used eq.(Eq). Note that it is set by the short- 
distance behaviour [x <C x c ) of the two-point correlation 
A^, since we actually have x = 0. On the other hand, 
for n^ < —3/2 the other terms T q with q > 2 are set by 
the long-distance behaviour of the correlation A^, i.e. by 
the scale R one is looking at by studying Sr. This is the 
case for the power-spectra P<p of cosmological interest, see 
eq.(43). Moreover, for large q the different factors A^ in 



eq.(|74|) are almost uncorrelated and it makes sense to use 
the approximation: 



q>2: T q ~\l [al(R L )Y 



(76) 



Besides, for a power-law power-spectrum P,p(k) as in 
eq. (|43|) we have (at scales larger than the coherence length 

X c )\ 



A^(xi,x 2 ) oc |xi -x 2 



_vJ-(»*+3) 



(77) 



Then, since n^ — —2.4 is close to —3 the correlation 
A0(xi,x 2 ) is almost constant over most of the integra- 
tion volume so that the approximation ([76]) should be 
quite reasonable. This yields for the trace T(Ao) defined 
in eq. (|7l|) : 

T(Ao) ^ 2iX {4> 2 )~2i\ al{R L )+\n [l + 2i\ al(R L )] .(78) 

Next, we could try to consider a saddle-point with re- 
spect to both <5_l(x) and A(x). However, one can check that 
this procedure does not give meaningful results. In fact, 
as described above in the derivation of the exact eq. ( |60| ) 
the generating function ipL(y) of the linear density field 
itself is not governed by a unique saddle-point. Indeed, we 
had to take into account the contributions from all real 
fields <5l(x) in the path- integral to get the Dirac func- 
tional which eventually provides the correct result. Then, 
in order to keep this structure while using the spherical 
saddle-points (]7(j) it is natural to try to approximate the 
path-integral (|5q) by: 



■0(y) 



1L,R L 

2tt 



dA e 2 



\ Tr ln(l+2iAA ) 



x e 



—ySn(5h,Rr)+i\o{(4> )+Sl,r l ) 



(79) 



Here Sr(Sl,r l ) = F{§L,R h ) is given by the spherical col- 
lapse solution of the equation of motions, while the matrix 
AA<£ is given by eq. (|73|) . Thus, we have replaced the inte- 
gration over the real fields <$l(x) and A(x) in eq. (|58|) by an 
integration over the real numbers 5l,r l and Ao. In other 
words, we have replaced the average over <5l(x) and A(x) 
by an average over spherical saddle-points of the form (|70[). 
In the case of the linear regime when <5(x) = <5l(x) we sim- 
ply have Sr — Sl,r l and Rl = R. Then we can perform 
the integration over 8l,r l which yields the Dirac function 
5d[Xq + iy}- The integration over Ao is now straightfor- 
ward and we recover exactly the rigorous result (pH) (this 
also justifies the normalization factor 2tt we introduced in 
eq.(|79|)). Note in addition that the structure of the cal- 
culation is also the same (the Dirac functional is simply 
replaced by a Dirac function). Thus, the expression (J79| ) is 
actually exact in the linear regime. We show in AppJAl how 
to rigorously derive the expression ( |79| ) from eq. (|5q) in the 
linear regime. We also explain how it can be derived for 
the fully non-linear density field in the quasi-linear limit. 
Thus, eq.(79) is actually exact to leading order in the limit 
a — > 0. More precisely, for the non- linear density field it 
yields the exact exponential cutoff of the pdf V{5r) which 
dominates the dependence on a but the normalization is 
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only correct up to a multiplicative factor of order unity 
which may depend on 8 R but not on a. This factor arises 
from the jacobians J a and J& in eq.( ATq). We shall come 
back to this point below in eq.(B3). 

As explained above, the factor Trln(l + 2iAA^) is not 
very convenient for practical calculations hence it is nat- 
ural to use the approximation ( f78| ) which was derived for 
such spherical fields A(x). This yields: 



^(y) 



dd L dX 



-yS R +i\8 L +i\o\ (R L )- 1 ln(l+2i\crl (R L )) 



(80) 



where we removed the subscripts of the integration vari- 
ables 5l and A. Note that in this approximation the term 
(eft 2 ) has disappeared. In the linear regime, where 8r = 8l 
and Rl = R, we can again perform the integrations over 
Sl and A which yields: 



$ 



,(y)~yal(R) 



I Hi 



2y<rl(R)) 



(81) 



Then, expanding the logarithm around y = we ob- 
tain the cumulants (8 q L R ) c from eq. (|62|) . This leads to 
{8l,r)c = and: 



q>2: (8 q ) c ^2^{q-l)\a 2q {R). 



(82) 



In particular, we obtain for the variance: u 2 {R) — 2cr4(_R), 



which must be compared with eq.(47) and eq. (J49J) . Of 
course, the deviation from the exact result (ES) is entirely 
due to the approximation (|76|). For a power-law power- 



spectrum P^ as in eq.(|43|) with n 
this deviation T>: 



4> 



-2.4 we obtain for 



V = 



2o A 



[/^^^( Xl ,x 2 )] ; 
where we used eq. (|77|) and 



1.11 



(83) 



the result (see 



Peebles fc Groth (1976)) 



dxi dx2 
V V 



Xi -x 2 



72 (2R)- 



(3- 7 )(4-7)(6- 7 ) 



(84) 



Thus, we see that the approximation ([76|) is quite reason- 
able. Note that it becomes exact in the limit n^ — > —3. 
Then, from eq.(82) we get for the skewness and the kur- 
tosis: 



D. 



(Sl,n 



(!■' 



2.83, D A 



(8 4 



L.R/c 

a 4 
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since we have: 



(85) 



(86) 



(Sl R ) c ^8a% (6 4 >R ) C ^A8 4. 

On the other hand, for the same power-spectrum (n^ = 
—2.4) Peebles (1999b) obtains by a direct numerical cal- 
culation (i.e. without the approximation (kq)): 



D 3 ~ 2.46, D 4 ~ 9.87 
and: 

(SI,r)c 



(87) 



8ctS 



1.02, 



(SU)c 



48 a5 



1.02 



Thus, we see that the approximation ( |76| ) is quite satis- 
factory. In particular, the behaviour of the high-order cu- 
mulants (8 q L R ) c is very well reproduced as soon as q <^ 3. 
Then, we can expect the expression (pQ) to provide a good 
approximation to ip{y), both for the linear and the quasi- 
linear regimes, since we do not add any approximation 
in order to describe the non-linear effects encoded in the 
relations <5l._r «-> 8r and Rl *-* R. 

Going back to eq. ( |80| ) , we can make the change of vari- 
able A — > —iz J o 2 ARl) which yields: 



4>(y) 



dS L f +icc dz e ~ ySrt+z[1+SL/rT i( RL ^ 



-oo °1(Rl) J-roo 27TJ y/T+2z 



(89) 



Then, after a change of variable and using the integral 
representation of the Euler Gamma function: 



(90) 



rto / 


+4 °° ^ _„ , 
- — z e 

-zoo 2m 




we obtain 








4>(v) - J 


-co 


d8 L 1 


1 + 


V^a 2 (R L ) \ 


X 


e~ 


y$R e -^[i+*i/ff? 


(Rl)] 



<(Rl) 



>0 



n-1/2 



*l(R L ) 



(91) 



where the factor 9 is the Heaviside function with obvious 
notations. Let us recall that in eq.(pTl) we have: Sr(Sl) = 
T{5 L ) andR L = R(l + 5 R ) 1 / 3 . 

4.4. ThepdfV L (5 L , R )andV(5 R ) 

In fact, since the generating function is written in an inte- 
gral form we can directly derive the pdf V(Sr) using the 
inverse Laplace transform (@) . Indeed, the integration over 
y simply gives the Dirac function 5d[5r — T{8l)\- After a 
trivial integration over S we obtain: 



V(S R ) 




1 



2wo%{R L ) F'(S l ) 

W+ s l/°1(Rl)} (92) 



where the quantity 6l is given by the condition: !F(8l) — 
8 R and R L = R^I + Sr) 1 ^ 3 . In particular, the pdf Vl(8l,r) 
of the linearly evolved density field is given by: 



Vl{8l,r) - 1 



Ql,r 
8l.r 



>o 



i 



-1/2 



2nal(R) 
e -J[H-*i.,fl/<^(fl)]_ (93) 



It is obtained from eq.(P2|) by setting J-{8l) — &l and 
Rl = R. For Gaussian initial conditions the pdf Vl,{8l ,r) 
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exhibits a specific scaling over the variable v = Sl, R /&(R), which yields: 
which contains all the time and scale dependence of the , 

pdf. As described above, within the approximation (um -p (§„) = Q \ \ 
we have a 2 ~ 2ci hence we get: V 

S L S l 



<{Rl) 



>0 



a(R) V2al(R)' 



(94) 



1 



Then, we see that the linear pdf (|93|) shows a scaling over 
v as in the Gaussian case: 



V L (S L , R )dS LiR = V^(u)du, 

but the pdf V l is now given by: 



vt\y) 



(l + V2i/>0) 1 



vT+75^ 



■i(l+V2t/) 



(95) 



(96) 



As noticed in paper I] , in the Gaussian case the non- linear 
high-density tail is well described by a simple spherical 
model. In fact, this model gives the exact exponential de- 
pendence of the pdf on the variance a, at leading order in 
the limit a — > 0, as was also shown in Valageas (1998) . Let 
us recall here how to build this simple model. It is based 
on the approximation: 



d<5 (1 + 5)V{6) 



dS L V L (S L ) 



(97) 



Sl.r l 



which merely states that the fraction of matter within 
spherical cells of radius R with a non-linear density con- 
trast larger than S R is approximately equal to the fraction 
of matter which was originally enclosed within spherical 
cells of radius Rl with a linear density contrast larger 
than 5l,Rl- Here Rl and Sl, Rl are related to the non- 
linear variables R and 5r by the usual relation (p5l), as 
in eq. (J92]) . Note that t his is similar to the usual Press- 
Schechter prescription ( Press fc Schcchtcr (1974)]), with- 
out the factor 2. Then, substituting the scaling (|95|) into 



eq.(97) and differentiating with respect to Sr we obtain: 



V s (6r) = 



>LM. L 



_J_ ^_ ?(")(„), with p = 

1 + S R dS R L V ; ' a(R L ) 



(98) 



Here the subscript "s" refers to the "spherical" model. In 
order to recover the variables used in the Gaussian case 
( paper Il| ) it is convenient to introduce the variable r given 
by: 



_ Ol,R 



a(R) 



<t{Rl) 



-ua{R), 



(99) 



which removes the dependence on the amplitude of the 
rms linear density fluctuation, and to define the function 
Q{t) by: 



G(t) = f[6 L , RL ] = S R 

which obeys again cq.(pi|). Then, the relation 



V S (S R ) 



1 + Sr a(R) \G'{t)\ 



Vt\v) 



(100) 
8|) writes: 

(101) 



\Q'{r) 



1 



<{Rl) 



-1/2 



2na 2 AR L ) 1 + Sr 



■h[i+SL/o$(RL)] 4 (102) 



We can see that we recover the expression (KJ3), except 
for the multiplicative factor where T'(Sl,r l ) has been re- 
placed by (1 + Sr)\Q'{t)\. Thus, as in the Gaussian case, 
we find that the simple spherical model ( p7J ) yields the 
exact exponential cr-dependent term of the pdf. However, 
the multiplicative factor may not be exact. Here, we note 
that the multiplicative factor which appears in eq.(|92|) 
is not exact either, as stated above, below eq. (|7^) , see 
also App.JAJ. In the Gaussian case, we c hecked by a com- 



(1998 



parison with numerical simulations (e.g., Valageas 
paper H| ) that the spherical model (J97J) provides good re- 
sults up to a ~ 1. We can expect a similar accuracy in the 
non-Gaussian case studied here. Therefore, the expression 
(102) should give a good description of the non-linear pdf. 
This means that the multiplicative factor which appears 
in cq.(102) gives a reasonable approximation of the term 
induced by the Jacobians which arise in the non-linear 
case, see eq.( A.15| ). 




Fig. 2. The pdf V{8 R ) for the isocurvature CDM sce- 
nario, with ri0 = —2.4 and a — 0.74. This corresponds to 
(70 = 0.70. The solid line shows the theoretical prediction 
from (102) for 7 , (<5h)- This corresponds to the "spheri- 
cal" model. The dotted line is given by eq.( 102 ) where 
no attempt is made to take into account the Jacobians 
which arise in the non-linear case. The solid line (spheri- 
cal model) should give better results. The dashed-line dis- 
plays the pdf Vl{Sl,r) of the linearly evolved density field, 
from eq.(p3|). The histogram shows the results of numeri- 



cal simulations by Robinson & Baker (2000) for the linear 
V<tiV L {5 L , R ). 

We show in FigJ| our results for n^ — —2.4 and 
a = 0.74. The dashed-line shows the pdf Vl(Sl R ) of 
the linearly evolved density field, from eq.(§)3). As could 
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be seen from eq.(93) it is strongly non-Gaussian, what- 
ever the value of a$. In particular, for large overdensitics 
$l,r the pdf exhibits a simple exponential cutoff (mul- 

1/9 

tiplied by a power 5 L ^ ) which is much smoother than 
a Gaussian cutoff. Hence the number of extreme events 
(i.e. Sl,r S> a) is much larger than for the Gaussian 
case, as can be checked by comparison with the lower 
dotted curve in Fig||. We can see that for Sl.r > 
our predictions agree very well with the results of nu- 
merical simulations from Robinson fe Baker (2000)| for the 
linear density field (shown by the histogram). Note that 
for the theoretical curve we used the exact value of a^ 
obtained from eq. (pa) in order to get a meaningful com- 
parison. Hence the variance of the approximate pdf ( |93| ) 
underestimates the exact result by a factor T> ~ 1.11. 
However, we can check that the agreement is quite good 
for 5lr > so that for pr actical purposes one can di rectly 
use eq.(P3[). For instance, Robinson fc Baker (2000) actu- 



ally us e the linear pdf to apply t he usual Press-Schechter 
recipe ( Press fc Schechtcr (1974) ) in order to estimate the 
mass function of just-collapsed objects. 

On the other hand, we note that our estimate (B3h fails 
for 8l,r Ss 0. In particular, we see that the pdf vanishes 
for 6L t R < —a?(R), that is the exact lower bound Sl,r > 
— {<t> 2 } has been replaced by 5l,r > —cr^(R). This is due 



to the approximation (|76|). Indeed, 
obtained from eq.(0) and eq.(pol) is: 



the exact linear pdf 



V L (S 



L,R) 



-\-ioc 



2iri 



2 ))-iTr ln(l+2»A v A*) 



(103) 



This expression clearly shows that the linear density con- 
trast obeys the lower bound Sl,r > — {<t> 2 )- Indeed, for 
&l,r < ~{4> 2 ) we can push the integration path over y to 
the right (i.e. Re(y) — » +oo) so that the integral vanishes 



(we have A y 



-iy). Note that this behaviour is rather 



different from the Gaussian case where the linear density 
contrast can take negative values of arbitrarily large am- 
plitude. In fact, the lower bound for 8l,r can be directly 
seen from eq.(J45|). Note that this bound is directly related 
to the time U at the end of inflation where eq.(Efi) holds. 



Indeed, from eq.rt44J) we see that 



D + (t) where 



the growing mode has been normalized by D + (ti) = 1. 
This minimum value of Sl.r, is also overestimated by 
the numerical simulations which only start at a redshift 
Zi ~ 39. Nevertheless, it is clear that the approximate 
bound 5l,r > — o\ gives the value below which under- 
densities become very rare (i.e. the pdf should decline for 
<>l,r ^ — crV). Besides, for practical purposes one is often 
mainly interested in the behaviour of overdensities <5r > 
so that the expression s (|92[ ) and (|93|) should remain useful. 
As noticed in eq.(Bq), the linear pdf given by eq.(B3|) 
exhibits a scaling over the one variable v defined in eq.(B4), 
as in the usual Gaussian case. This scaling property was 
actually used in order to obtain Pl(Sl,r) from the numer- 
ical simulations. However, this property is not exact and 
it fails for underdensities since as noticed above the linear 
density contrast also satisfies the lower bound — (</> 2 ) which 



does not scale as a (e.g., contrary to a it is independent 
of scale) . 

Next, we show the theoretical prediction V(Sr) from 
eq.(j92|) (dotted line) and from eq.(102) (solid line) from 
the spherical model ( p7J ) for the actual non-linear density 
field. In Fig.H we used the approximation: 



H&) 



5 
L5 



(104) 



which has been shown to provide a very good fit to the 
exact spherical collapse solution for all values of f2 m and 
f^A, see Fig. 2 in Bcrnardcau (1994b) . This means that the 
dependence on O m and ^a is negligible (in the quasi-linear 
regime where we do not consider virialized objects) and 
the pdf only depends on the rms fluctuation a (and the 
slope riff, of the power-spectrum) . Unfortunately, there are 
no available results from numerical simulations to com- 
pare with our prediction. However, in view of the reason- 
able accuracy of our prediction for the linear pdf Vl(5l,r) 
and the good results obtained for the Gaussian case (see 
paper II) we can expect the expression ( |l02j ) to provide a 
good estimate of the exact pdf V(Sr). It should perform 
better than eq.(J92j) where no attempt was made to take 
into account the Jacobians which arise in the non-linear 
case, see App.|A|. Of course, we recover the usual features of 
the non-linear evolution which increases the high-density 
tail of the pdf. Note again the much larger number of high 
density events with respect to the Gaussian case with the 
same variance (see Fig.pl). 

5. Standard slow-roll inflation 

5.1. Initial conditions 

The two classes of models we investigated in Sect.|| and 
Sect.[| are strongly non-Gaussian (except for a ~ 2 in 
the first case) and they apply to "non-standard" scenar- 
ios for the generation of primordial density fluctuations. 
For instance, the second model is based on a multifield 
inflationary scenario. By contrast, in standard slow-roll 
inflation the primordial density perturbations should be 
very close to Gaussian. Nevertheless, they may still show 
some small deviations from Gaussianity. In particular, one 
is led to consider models where the perturbed primordial 
gravitational potential ^(x) is given by: 



fc(x) = 0(x) + e (0(x) 2 - (4> 2 )) 



(105) 



where 4>(x) is a Gaussian random field with zero mean and 
power-spectrum P^,(k), as in eq.(|42f). Such models arise in 
standard slow-roll inflation if we keep track of the per- 
turbations of the inflaton up to the second-order (e.g., 
pangui et al. (1994)| , [Falk et al. (1993| . In fact, eq.([T05|) 
can be seen as the first two terms of a Taylor expansion 
so that it applies to a large class of slightly non-Gaussian 
models. Thus, the parameter e is taken to be small so that 
the primordial density fluctuations are close to Gaussian. 
Note that e has the dimensions of 6~ l . We shall define 



below what we mean by "e being small" (see eq.( |l38| )) 



P. Valageas: The quasi-linear regime for some non-Gaussian initial conditions. 



13 



Then, our goal here is to obtain the pdf V(5r) of the 
density fluctuations in the quasi-linear regime up to first 
order in e (to go beyond this order we should first take 
into account the possible higher-order terms in eq.(105)). 
Note that for e = we must recover the case of Gaussian 
primordial density fluctuations studied in paper II. The 



gravitational potential $(x) is related to the primordial 
linear density fluctuations Sl (x) by the Poisson equation 
(in comoving coordinates): 



A$ 



^-■KQp 



<Ux) 



where a(t) is the scale-factor and ~p the mean comoving 
density. Taking the Fourier transform of eq.(106) one can 
eventually write the linear density field 8l (x) at the time 
of interest as: 

6 L (x) = J dx' W(x, x') [0(x') + e (<Xx') 2 - (cf> 2 ))} (107) 

where the kernel W(x, x') is given by: 

W{x,y>) = -\^H-*D + (t) J '^ e * k -(— VT(fc) 

(108) 

Here we introduced the Hubble constant Hq and the den- 
sity parameter fi m today, while D+(t) is the usual lin- 
ear growing mode at the time of interest, normalized by 
D + (to) = 1 today. The function T(k) is simply the adi- 
abatic CDM transfer function (normalized to unity for 
k — > 0). Indeed, contrary to the isocurvature model de- 
scribed in Sect.[| we need to take into account the de- 
viations of T(k) from unity on the scales of interest for 
large-scale structure formation. The factor k 2 arises from 



the Laplacian in the l.h.s. in eq.(106). Note that the kernel 
W(yi, x') is homogeneous, isotropic and symmetric since it 
only depends on |x — x'|. Hence for any real fields f\ and 
/ 2 we have: 



(109) 

Thus, eq. ( |l07|) defines the initial conditions of our system. 
We can note that the pdf Vl(3l,r) of the lin- 
early evolved density field was already studied in 
Matarrese et al. (2000) using path-integral methods like 



those described in Sect. 5.2 below. However, they did not 
investigate t he effects of the non -linear dynamics. On the 
other hand, Verde et al. (2000) studied the observational 
tests which may constrain such deviations from Gaussian 
initial conditions. 

5.2. Generating function 

In order to derive the pdf V(Sr) we again introduce the 
Laplace transform ip(y) as in eq.(||), which yields again 
eq.(52). Next, we introduce the linear density field through 
the Dirac functional 6d[Sl(x) — J dx' W(x, x'){0(x') + 



e(</>(x') 2 — (</> 2 ))}]. This Dirac functional can again be ex- 
pressed through an auxiliary real field A(x) so that the 
analog of eq.(154) now reads: 



fl>(y) = y[d0(x)][d ( 5 L (x)][dA(x)]e-^ [5il -^- A * 1 ^ 

ol A.(<5 i -W.[0+e(0 2 -<0 2 ))]) ( U0 ) 



x e 



where we defined the inverse A . of the two-point corre- 
lation A^ of the Gaussian random field <p, which is again 



(106) given by eq.(|46|). The mean (4> 2 ) is again given by eq. 



and it does not depend on the field </>(x) over which we 
integrate. As in Sect.|l|we do not write explicitly the nor- 
malization factor of the path-integrals. Next, we introduce 
the kernel Aw(xi,x 2 ) defined by: 



A w (xi , x 2 ) = S D (xi - x 2 ) / dx A(x) W(x, xi ) 



(111) 



so that eq.(llC) reads: 

f(,(y) = /[d0(x)][d,5 L (x)][dA(x)]e-^ [ ' 5il -^- A * 1 



x e 



i\.S L -i\.W.cf>-ie<p.A w .<f>+ie(4> 2 )(\.W.l) 



(112) 



As explained in Sect. 5.1 we only consider the first-order 
term in e. Therefore, we expand the exponential in 
cq.( |ll2| ) up to first-order in e. This yields: 

My) = /[d^(x)][d ( 5 L (x)][dA(x)]e- y5 ^ a - W/ ^-^' A * 1 ^ 



x e iX - 5L [l + ie{<j) 2 ){\.W.l)-ie{cj)A w .(j))\ .(113) 

This path-integral is Gaussian over the random field (j). 
Thus, we first mak e th e change of variable = 0' — 
iA<p.W.\ so that eq.(113) writes: 



ip(y) 



[d0'(x)][d<5 L (x)][dA(x)] e 



-y6 H +iX.SL-iX-W.A^.W.X 



ih-W.h) = (/2.W./1) = / dx 1 dx 2 M/(x 1 ,x 2 )/ 1 (; Cl )/ 2 ( ;C2 ). x e 



■AI 1 .*' 



1 + ie(cf> 2 ) (X.W.I) - ie{<j>'A w 4') 



-2e{4>'A w .^.W.\) + ie(\.W.A (t ,.A w .A l p.W.\) 



(114) 



Here we used the fact that all kernels W, A^ and A\y are 
symmetric. Then, using Wick's theorem we can perform 
the Gaussian integration over </>'(x) which yields: 



i>(y) 



<Wi(x)][dA(x)] e-vSn+iX-SL-^-w.At.w.x 

1 + ieiX.W.A^.Aw.A^.W.X) 



(115) 



where we used A<ji(x, x) = ((f) 2 ). Next, we define the two- 
point correlation Aio(xi,x 2 ) of the linear density field in 
the case e = by: 



A L0 (xi,x 2 ) = (<5 l (xi)(5l(x 2 ))o 



((W.0)(xi)(W.^)(x 2 )) 

(116) 
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where we used eg. (107 ). H ere the subscript "0" refers to 
"e = 0". The relation (116) also writes: 

A L0 = W.As.W 



(117) 



The path-integral (115) is Gaussian over the field A(x). 
Thus, making the change of variable A = A' + iA^g .Sl and 
integrating over A'(x) using Wick's theorem we finally get: 



(118) 



1p(y) = / [<Wi(x)] e -vSnlS L ]-i6 L .A L0 .6 L 

l + e(A7 1 .Ty- 1 .5 L ).(W" 1 .5 L ) 2 



where we used /dx 8l(x) = 0. Here we introduced the 
short-hand notation (W~ 1 .8l) 2 for the vector: 



(W-KS L ) a (pn) = (7dx' T^- 1 (x,x')<5 L (x'; 



(119) 



Note that for e = we recover the expression ([14]) derived 
for Gaussian initial conditions. This was to be expected 
since for e — the linear density field is actually Gaussian. 
Indeed, in this case we have 8l =W.<j> and </> is a Gaussian 
random field. Of course, the procedure we described above 
can be extended up to any order in e (provided we know 
the initial conditions up to the required order). Indeed, by 
expanding the exponential which appears in eq.( |llO| ) up 
to the needed order in e the integrations over the fields 
</>(x) and A(x) are Gaussian and can be easily performed 
using Wick's theorem. Then, one eventually obtains an 
expression of the form (118), where the term in brackets 
is an expansion over e up to the required order. Note that 
this is similar to the standard Edgeworth expansion where 
one directly expands the Laplace transform ip(y) around 
the Gaussian value ipo(y) — e y a ^ (see any textbook on 
probability theory). Therefore, the expansion over e does 
not necessarily converge: it may only provide an asymp- 
totic series (e.g., pramer (1946)). 



5.3. Steepest-descent method 



We now need to evaluate the pa th-integr al (118). As 
for the Gaussian case studied in paper II) or the non- 
Gaussian generalization described in Sect.fil we can use 
a steepest-descent method for the quasi-linear regime. We 
refer the reader to paper II for a detailed presentation of 
this method. First, we define the rescaled generating func- 
tion i/;(y) by: 

^(y) = V (y/a (R) 2 ) (120) 

where cr (R) is the variance of the linear density field for 
e = 0: 

dx! dx2 



MR) 2 = Kr)o = 



V V 



Al (xi,x 2 ) 



(121) 



Using eq.(118) we get: 



i>(y) 



[«Wi(x)] e 



-S[S L ]/a (Rf 



1 + e (A7 1 . ^.SlUW- 1 ^!.)' 



where we introduced the "action" S[5l] given by: 
S[S L ] = y Sr[6l] H -z — Sl-A-lI-Sl 



(123) 



In the quasi- linear limit <7o — *■ the path-integral (122) is 
governed by the minimum of the action S. As shown in 
paper II this spherically symmetric saddle-point is given 
by: 



<5l(x) = 8 LtR 



dx' A L0 (x,x') 



^ m. m,A T ,n (0{x'<R L )\ 



MR 



(W.A^.W). 



V L 



(124) 



where x' is a dummy variable. The variable 5l,r l IS given 
by the implicit equation: 



7L,R L 



-y 



f'[8 LtRL ]a 2 (R L )/a 2 (R) 



1-^[8l,r l ]4tS l , r 



1 dcrp 

iL <J (Rl) dR 



mm 



.(125) 



Here the function J 7 (Si Jt ii L ) is the usual spherical col- 
lapse solution while Rl is the Lagrangian mass scale, see 



eq.(|25|). These results were also used in Sect. 3. 2 where 
we investigated a simple generalization of the Gaussian 
weight. Next, we introduce the functions t(Sl) and Q(t) 
by: 



-(Sl) = 



- S L a (R) 



and: 
G{t)=F[8 l {t)) = 5 r 



(126) 



(127) 



as in eq.(E9J) and eq.(J30j). Then, we define the function 
<po(y) = toxd.S[8l] as the value of the action S at this 
spherically symmetric saddle-point (i.e. the minimum of 
the action) . As shown in paper II it is given by the implicit 
system: 



■V0'(t) 



<po(y) =yQ{r) + y 



(128) 



in a fashion similar to eq.(32). T hen, at leading order in 



do we can write the path- integral ( 122 ) in the quasi- linear 
limit as: 



i/j(y) = e ~<Po(y)/MR) 2 



1 



MRl 



a (Rr a Q (R L ) 3 



where we defined: 



J 3 {R L ) = W. 



(X! < R L ) 



A A .W. 



0(x 2 < R L ) 
Vl 



(129) 



.(130) 



Here xi and x 2 are dummy variables. The normalization 
of ip(y) in eq.(129) is set by the constraint i/j(0) = (1) = 1. 
Next, we define the dimensionless quantity e by: 



(122) ~<Rl 



MRl) 
MRl) 3 



(131) 
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It is convenient to express the quantities Co and J 3 in 
terms of the power-spectrum P</,(k) of the Gaussian ran- 
dom field 4>(x). This yields: 



<j (Rl) 2 = (2tt) 6 fdkP^k)W(k) 2 F(kR L ) 2 



(132) 



and: 



J 3 (R L ) = (2ir) 9 / dkdk' P^{k)P^k')W{k)W{k') 

x W(|k + \<!\)F{kR L )F(k'R L )F{\k + k'\R L ). (133) 

Here W(k) = W^dkl) is the Fouri er tr ansform of the win- 
dow function W(x) defined in eq.(108) (i.e. it is related to 
the kernel W by W(xi,x 2 ) = W(x 2 - x x )) while F(fc.R) 
is the Fourier transform of the top-hat of radius R: 



F(kR) 



V 



o sm(kR)-(kR)cos(kR) 

6 (kRj 3 ' { ' 



Note that in eq.(134) we did not introduce the factor 
(27r)~ 3 used in the definition (||) of the Fourier transform, 
in order to obtain the usual top-hat window F(kR). As in 
the previous sections we approximate the power-spectrum 
Po(k) of the linear density fluctuations (when e = 0) by 
a power-law. Since we have P<f>(k) oc k~ 3 (i.e. a scale- 
invariant primordial power-spectrum, with small and large 
scale cutoffs) for the standard inflationary scenarios we in- 
vestigate here this means that W(k) 2 oc k n+3 on the scales 
of interest. Indeed, the power-spectrum Po(k) is related to 
P*(k) by: 



P (fe) = (27rfW(k) 2 P 4> (k), 



as can be seen for instance from eq.( |l32j ). Then, from 
eq.(132) and eq.(|l33j) we get: 



a (R L ) 2 <x K 



-} n+3 \ J 3 (i? L ) oc i?- 3( " +3)/2 



(136) 



which implies that the quantity £(Rl) defined in eq.([13lj) 
does not depend on Rl, within this approximation. 



Moreover, from eq.(132) and eq.(133) we sec that we have 
the order of magnitude estimate: 



dkP (fc) 



1/2 



vw 



(137) 



Going back to the definition ( |105|) of the initial conditions 
we see that we can write the primordial perturbations of 
the gravitational potential as: 



$(x) 



(x)+0(l)e 



(x) 2 



vw 



(138) 



where 0(1) stands for a numerical factor of order unity. 
It is clear from this expression that for a scale-invariant 
primordial power-spectrum the initial conditions are close 
to Gaussian if e <C 1. Hence e is the relevant dimensionless 
parameter which describes the amplitude of the deviations 
from Gaussianity and the usual inflationary scenarios lead 
to e< 1. 



Here it is interesting to consider the linear regime, that 
is the statistics of the linear density field. Then, we have 
Q{t) = -t and r = y = -5 l ,r = Sr, so that ip (y) = 
-y 2 /2. This yields from eq.fll29|): 



i>L{y) 



aV^l/Z 



[l-f«] 



(139) 



where the subscript "L" refers to the "linearly evolved" 
density field. Then, from the expansion (0) we obtain for 
the third-order moment of the linear density field (<5£ R ) 
and for the skewness D3: 



(S 



l,r) 



'0) 



(4 



= 6eal £> 3 = 



(Sin) 

(SI,r) 3/2 



6e (140) 



at first-order in e. Thus, we see that the parameter e is 
directly related to the skewness of the linear density field 
on the scale of interest, due to the slight non-Gaussianity 
of the primordial density field. 

5.4. The pdf of the density contrast 

From the rescaled Laplace transform ip(y) obtained in 
eq.(129) we can derive the pdf V(Sr), using the inverse 
Laplace transform (H). This yields: 



V{8 R ) = 



+ ZOO 



dy 



2TTi(T 2 {R) 



,[y6 R -<po(y)]/<J% 



r 31 



(J, 



(141) 



where the function r(y) is given by eq.(128). The pdf 
T > l{5l,r) of the linearly evolved density field is obtained 
from eq.(141) by using <fio(y) = —y 2 /2 and r(y) — y which 



(135) gives: 



V L (S 



L.R) 



-\-ioc 



dy 

2ni(jQ 



olvS l, R+y 2 / '2] jcl 



e —5 



.(142) 



Then we can perform the Gaussian integration over y in 
eq. (|l42j ) which yields (see |Gradshteyn fc Ryzhik (1965) , 
§3.462.4): 



Vl{5l,r) 



-&\ tR l{iol) 



2na (R) 



l + i2-^ 2 H 3 f-£M 



V2, 



(TQ 



(143) 



where H3 (x) is the Hermite polynomial of order 3 defined 
by: 



H 3 {x) = 8x 3 - 12a;. 



(144) 



Finally, we display in FigJ3 our results for slightly non- 
Gaussian initial conditions: e = 0.15. In order to compare 
the pdf with the case of Gaussian initial conditions we also 
show the results obtained for e = (dotted curves). Let 
us recall that eq. fll42j ) for the pdf of the linearly evolved 
density field is exact to first-order in e. In agreement with 
cq.( |140| ) we can check that e > leads to a small pos- 
itive skewness for the linear density field. In particular, 
the high-density tail of the pdf is slightly enhanced. We 
can see that this feature remains valid for the non-linear 
pdf V(Sr). However, as gravitational clustering proceeds 
the deviations from Gaussianity become dominated by 
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Fig. 3. The pdf V{5r) for n = -1 and cr = 0.74, 
for slightly non-Gaussian initial conditions: e = 0.15. 



The solid line shows the prediction of eq.(141) for the 
non- linear pdf V^r) while the dashed-curve shows the 
pdf Vl (Sl,b) °f the linearly evolved density field from 
cq.(143). The dotted lines display the results obtained 
with e = 0, i.e. for Gaussian initial conditions. 

the non-linear dynamics over the slight primordial non- 
Gaussianity. More precisely, we can obtain the low-order 
moments of the density field through the expansion (Q). 
At the lowest order in y we have r = y + 0(y 2 ) and 
<p Q (y) = -y 2 /2! + S 3 y 3 /3\ + 0{y 4 ) where the parameter S 3 



is given by the standard result (e.g., Bcrnardcau (1994a)| ): 
34 

y 



s, = 



(n + 3). 



(145) 



Then, we obtain from eq.(129) the expansion 

1 



i>{y) 



i 



2 2 

-,y CT o 



i 



3„4 



r S 3 y% 4 



«/ 3 <7o 



2\ a u 3! 
This yields for the lowest order moments: 



0(y 4 



K)=4> (S R )=S 3 a 4 +6 



? ~3 



and for the skewness D 3 



D <*«> -6£ 



S3 c . 



(146) 



(147) 



(148) 



This relation clearly shows that the relative importance 
of the primordial deviation from Gaussianity is greater at 
earlier times before gravitational clustering builds up, as 
could be expected. Note that in Fig.|| we used the gen- 
erating function <po(y) defined by the branch of eq.(128) 
which is regular at the origin y — 0. Since here we con- 
sidered the case n — — 1 this function actually shows a 
branch cut along the real negative axis for y < y s (with 
y s < 0). Then, for large density contrasts one needs to 
take into account the second branch of (fo(y) defined over 
the range y s < y < which is singul ar at the origin. 
However, as shown in Fig. 6 in paper II this is irrelevant 
for the c ase displa yed in Fig.|| for Sr < 7. We refer the 
reader to |paper II for a detailed discussion of this point. 
Of course, for the non-Gaussian model we study here one 



simply needs to follow the procedure outlined in paper II 
for Gaussian initial conditions and to add to the relevant 



expressions the term in the brackets in eq.(141) which de- 
scribes the effects due to the small primordial deviation 
from Gaussianity (to first order in e). 



6. Conclusion 

Thus, using a non-pcrt urbative method developed in a 
previous work ( paper II ) , we have described in details how 
to obtain the pdf P(6r) of the density contrast within 
spherical cells in the quasi-linear regime for three specific 
non-Gaussian models. 

The first case is a straightforward generalization of the 
Gaussian scenario and it can be seen as a phenomenolog- 
ical description of a density field where the tails of the 
linear pdf are of the form V(5l) ~ e - '" 51 -' , where a is 
no longer required to be equal to two (as in the Gaussian 
case) . Then, we have shown that the derivation of the pdf 
presented for the Gaussian case can be directly extended 
to this model. This provides again exact results for the 
pdf in the quasi-linear limit. 

The second scenario is a physically motivated 
model of isocurvature cold dark matter presented in 
It arises from an inflationary scenario 



Peebles (1999a 



with three scalar fields. This case is slightly more diffi- 
cult as one needs to adapt the method to this specific 
model. Moreover, in order to get simple analytical re- 
sults one must introduce a simple approximation (which is 
not related to the gravitational dynamics but to the non- 
Gaussian properties of the initial conditions). However, 
we have shown that even with this approximation we get 
good results for the linear pdf Pl(Sl,r) for <5l,_r ^ by 
comparison with numerical simulations. For lack of data 
we could not check our prediction for the non-linear pdf 
V(Sr) but we can expect a good agreement of the same 
accuracy. 

Finally, the third scenario corresponds to the small pri- 
mordial deviations from Gaussianity which arise in stan- 
dard slow-roll inflation. We obtained exact results for the 
pdf of the density field in the quasi- linear limit, to first- 
order over the primordial deviations from Gaussianity. 

This study shows that our approach is powerful enough 
to be applied to a large variety of initial conditions. Note 
that our predictions for the linear pdf Vl(Sl,r) for such 
non-Gaussian scenarios may be used to estimate the mass 
function of just-collapsed objects, using a straightforward 
extension of the Press-Schechter recipe. Note that we shall 
discuss the Press-Schechter prescription in the light of the 
formalism developed in paper II and used in this work in 



a companion article (paper IV). 

To conclude, we note that our approach being non- 
perturbative it can in principle be applied to the non- 
linear regime. Indeed, it does not rely either on the hy- 
drodynamical description. We shall prese nt a study of 
this non- linear regime in a future work, see paper IV , for 
Gaussian initial conditions. However, it is clear that this 
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can be extended to the non-Gaussian models described 
here. 

Appendix A: Isocurvature scenario. Reduction to 
ordinary integrals 



where the subscript "L" refers to the "linearly evolved" 
density field (so that 8r — Sl.r — &a) and we did not 
write the normalization constant of the integrals. Then, 
the integration over the variables a q with q > 1 yields the 
Dirac functions 8o{b q V). The integration over b r is now 
straightforward for r > 1 and we obtain: 



In this appendix we show how to derive the expression ( |79[ ) 
from eq.(pq). We first consider the linear regime where 



^l{v) 



eq.(79) is exact. The path-integral ( J58| ) involves the func- 
tional measure [d<5t(x)]. The latter can be defined from 
a discretization procedure (i.e. a spatial grid for the co- 
ordinate x with an infinitesimal spacing) but an equiv- 
alent formulation (e.g., Zinn- Justin (1989)) is to expand 
the function 8l (x) on a complete set of real orthonormal 
functions f q (in the Hilbert space C 2 ): 



da db e -y-o+iVaob +i^ )b v 

) 

x e -iTrln(l+2iAA^) 



(A.; 



*L( X ) =^2 a q /<?( X )- 
9=0 

The vectors f q obey the relation: 
/ y / g (x)/ r (x) = 5 q . r 



(A.l) 



(A.2) 



where the matrix A is obtained from A(x) = b 9(x < 
R). Then, using ao = 8l,r from eq.( |A.5| ) and making the 
change of variable \q = b^V we obtain eq.(f79|) in the linear 
regime (i.e. 8 l . Rl = S l ,r = 5 R and R L = R). 

Now, we investigate how the previous derivation is 
modified when we study the non- linear density field <5(x). 
We can again d efine the functional measure [d^(x)] (and 
[dA(x)]) by eq.(A3). However, we now introduce the new 
variable Sl,r l , which also defines 5r = T(8l,r l ) and 
Rl = R{1 + Sr) 1 / 3 . Then, we define the functions f q (x) 
by: 



which expresses the fact that they form an orthonormal 
basis. Here 8„ r is the usual Kronecker symbol. We intro- 



/;w = V^ /9 (^ x 



(A.9) 



duced a factor 1/V in the scalar product ( A.2 ) so that 
the functions f q are dimensionless, but this is not essen- 
tial (here V is a constant). The n, the functional me asure 
[d(5/,(x)] can be defined as (e.g., Zinn- Justin (1989) ): 



so that we still have the orthonormalization property: 



J q,r • 



[dS L { X )]=Ul[da q 



(A.3) 



Then, we write the linear density field <$l(x) as: 



g =o 



where Af is a normalization constant. In order to get a 
discrete basis in eq.( |A.l[ ) we restricted the density field 
to a large finite volume V. However, this is not essential: 
we choose V to be much larger than any relevant length 
scale and we can eventually take the limit V — > oo. Next, 
we can choose for the first basis vector /o the top-hat of 
radius R: 



,/ (x) = 9(x < R) 



(A.4) 



where 0{x < R) is th e us ual top- hat w ith obvious nota- 
tions. Then, from eq.(Al) and eq. ([A.2| ) we have: 



dx 

"ii = / TT /o(x)M x ) = h,R- 



(A.5) 



We can expand the auxiliary field A(x) on the same basis: 

CX) 

A(x)=^6 g / g (x) (A.6) 

q=0 



so that the path-integral (58) now writes: 



Mx) = 5>;/;(x) 

which also implies: 
' fib X 

y ~v lrl ' 

This yields: 

oo 

8l(x) = Sl,r l 0(x < Rl) + Y, < / 9 ( x ) 

9=1 

while we write the auxiliary field A(x) as: 

oo 

A(x)= £X /£(*). 

9=0 

Then, the path-integral (pq ) now writes: 

/oo oo oo 
d8 L M L l[da q l[db' r J a J b 
-OO i n 



(A.10) 



(All) 



(A.12) 



(A.13) 



(A.14) 



1-ySn 



^l{v) 



\\ da q JT db r e- yaa+l 

30 9=0 r=0 

i(0 2 )(A.l)-iTrln(l+2iAA^,) 



V J27=o a " b " 



x e 



9=1 r=0 



(A.15) 



x e 



(A.7) 



where we made the changes of variables {a q } q >p — v 
{8L,R L A a ' q }q>i} and{6 fl } > -> {b' q } q > . Ineq. (|A.15|) the 
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jacobians J a and Jb of these transformations are given by 
the determinants: 



Ja 



Det 



da,. 



OS 



LMl 



dai 

da'. 



, Jb 



Det 



db'. 



(A.16) 



These jacobians do not depend on the amplitude a(R) of 
the density fluctuations. On the other hand, the quasi- 
linear regime corresponds to the limit a — > for a fixed 
finite y (which leads to finite 8l,r l and <5r, as shown in the 
Gaussian case in paper II). In this limit, the generating 
function ip(y) is governed by the exponent in eq.(A.15) 
which depends on a and sets the cutoffs of order a for the 
density contrast <5r. Then, we can neglect the Jacobians in 



eq.(A.15): they do not contribute to ln(-0) at leading order 
for a — ► (however, they give a cr-independent prefactor 
which appears in ip, see the discussion below). Moreover, 
if we make the approximation <5r[<5l] — J : {5l,r. l ) we can 
again integrate over a' for q' > 1, which gives the Dirac 
functions 8d(P'V), and we eventually recover eq.(|79|). 

Let us now discuss the approximation <5r[<5l] ~ 
^F{5l,r l ) used above. The point is that in the limit a — > 



the path-integral ( p8| ) and the ordinary integrals ( A. 15 ) 
are dominated by the saddle-points of the exponent. As 
discussed in Sect. 4. 3, because of the spherical symmetry 
of the physics we investigate there exist some spherical 
saddle-points, of the form (J70|) . Then, if there are no other 
saddle-points (or if some other saddle-points exist but they 
yield a lower contribution to integral, which will vanish in 
the limit a —> 0) the path-integral is given at leading-order 
by its value for these spherical saddle-points. Note that the 



case of Gaussian initial conditions studied in paper II ex 



hibits the same behaviour in a simpler manner. Indeed, in 
this case there exists only one spherical saddle-point for 
both the linear and non-linear generating functions ip^ (y) 
and ip(y). Then, in order to keep only the leading order for 
a — ► one defines the rescaled generating function ip(y) 
by: 



i>(y) 



— e ~v{y° )l" 



(A.17) 



as in eq.(^9[), which can be expressed through the path- 
integral: 



-<p(v)/v — 



[cWl(x)] e 



-S[S L ]/a 2 



(A.18) 



in a fashion similar to eq.(EG). Next, in the quasi-linear 
limit a — > the path-integral in the r.h.s. in eq.(A.lE) is 
dominated by one saddle-point and taking the logarithm 
of both sides in eq.(A.lS) one obtains in this limit: 



difference is that in this non-Gaussian case there exists an 
infinite number of spherical saddle-points with respect to 
Sl(x). Therefore, at leading order the generating function 
ip{y) is not given by a unique contribution but rather by 
the sum of all contributions due to these saddle-points 
which are parameterized by the variables &l,r l and Ao- 

Finally, we note that the minimum of the action S[(5x,], 
or the value of the exponent in eq.(A.15), only pro- 
vides the leading order contribution to ln(^>) in the limit 
a — » 0. This means that it yields the exact exponential 
cr-dependent term of the pdf V(6r) which describes the 
cutoffs of the pdf for density contrasts which are large 
relative to a. In our case, this is the factor e~' 1+ ° L ' a '<t>>' 2 
which appears in eq.(62). However, it does not give the ex- 
act multiplicative prefactor of the generating function ip 
an d of th e pdf V(Sr). Indeed, the Jacobians J a and Jb in 
eq. (|A.1E| ) (or the Gaussian integration around the saddle- 
point in the Gaussian case) contribute to this term. In the 
case we study here, this means that eq. (|79f ) differs from the 
exact result by a multiplicative factor which may depend 
on 6l and A. Therefore, the expression (92) obtained for 
the pdf is correct up to a multiplicative factor which may 
depend on Sr but not on a. This point is also discussed in 
Sect. 4^4 where we compare eq.(£^) with the result (102) 
of a simple spherical model. The latter can be seen as an 
attempt to model this prefactor. Note that for the pdf of 
the linear density field we do not encounter this difficulty 
since, as shown above, eq.( [A.8| ) is actually exact. 
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